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Abstract 

In this article, we describe a recursive method to construct a subset of relevant Slater- 
determinants for the use in many-body diagonalization schemes that will be employed in our 
forthcoming papers on the simulation of excited many-body states in discrete electronic nanosys- 
tems. The algorithm is intended for the realistic simulation of nanodevices which typically requires 
the consideration of a large number of single-particle basis states (typ. 256) and the calculation of 
a sufficient number (typ. a few 1000) of relevant excited many-body states. 
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I. INTRODUCTION 



Realistic microscopic modeling of application-relevant electronic nanostructures typically 
requires the consideration of 100-1000 single-particle basis states (orbitals). Normally, due to 
the lack symmetries, such systems are inhomogeneous and anisotropic and furthermore may 
exhibit a long-range interaction. With a suitable unitary transformation and a subsequent 
division of the single-particle Hilbert space, the number of relevant single-particle states that 
require a many-body treatment of the Coulomb interaction typically can be reduced to <256. 
Even with this reduction, the resulting Fock space dimension obviously grows beyond prac- 
tical numerical limits, at least for a full exact diagonalization (full-CI). In order to address 
this problem one can take the path to consider a highly idealized model (such as Anderson- 
or Hubbard-like models , mapping the actual nanosystem to a reduced (albeit non- 

trivial) model system with a small number of effective parameters. Such an approach can be 
justified in the sense that only a small number of degrees of freedom typically are relevant 
for the investigated fundamental effects. Nevertheless, these relevant degrees of freedom and 
their associated effective parameters in general are non-linear functions of the actual exper- 
imentally accessible parameters (such as gate voltages and general external electromagnetic 
fields) and not known a-priori. As an advantage however, due to the simplified structure of 
the model Hamiltonians, one is able to understand the physics of such idealized models in 
detail (in particular, the possibility to correctly include continuous semi-infinite contacts). 
However, in order to determine effective parameters (as a function of external quantities 
and material properties) and to obtain a more realistic description of the nanosystem (in 
particular for the simulation of application- relevant structures) , methods which incorporate 
more microscopic details are mandatory. While ab-initio methods in principle provide this 
kind of information, they often involve approximation schemes (such as mean-field approxi- 
mations) that are not suitable for the description of many-body correlations which certainly 
cannot be neglected in nanosystems with quantum effects. Nevertheless, ab-initio basis 
functions and matrix elements might serve as input parameters for subsequent many-body 
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librium, idealized model interaction or periodic systems, limited-particle excitations, small 
single-particle basis, restriction to systems with special symmetries (such as spin and total 
angular momentum, leading to block- diagonal Hamiltonians) , restriction to Fock-Darwin 
single-particle basis (or even 2D-harmonic potential only). 

In the following, we discuss a recursive scheme which is able to describe a finite (or quasi- 
isolated) nanosystem with a large number of single-particle basis states (typ. 256) and to cal- 
culate a sufficient number (typ. a few 1000) of relevant excited many-body states. Note that 
the states we are interested in are not necessarily close to the groundstate. This approach 
consists of two steps: (i) Recursive construction of a large set of relevant Slater-determinants 
("bucket brigade algorithm") as discussed in this paper, (ii) Many-body calculations within 
the resulting Fock-subspace. 

II. RECURSIVE CONSTRUCTION OF RELEVANT SLATER-DETERMINANTS 

For an arbitrarily given single-particle orthonormal (ON) basis B = {4>j} the exact many- 
body Hamiltonian of the finite (i.e., discrete) system in second quantization reads as 

oo oo oo k—1 oo I— I 

H = Y.Y1 ^jkc]ck + yjkimclclcicm, (1) 

j=0 k=0 k=l j=0 1=1 m=0 

where ejk and Vjkim denote single-particle and Coulomb two-particle matrix elements of the 
total system Hamiltonian, and Cj is the annihilation operator for single-particle state j. 

Let Nmax be the maximum number of single-particle states that shall be considered. 
Hence, we consider a subspace spanned by B]^^^^ = {(pj | j = 0, Nmax — I}- The recursion 
procedure can be formulated in terms of a sequence of sets Sj^n of Slater-determinants. 
Here, J = 0, Nmax denotes the recursion step number and = 0, Nmax is the particle 
number. Sj^^ is defined via the following "bucket brigade" recursion as visualized in Fig. 

• Start (initial conditions) 

For J = 0,... , Nmax-. Sj,o = {\0)}, 

ioi N =1,... , Nmax-. So,N = ^- 

• Recursion step J —>■ J + 1 

For J = 0, Nmax - 1 and iV = 1, Nmax ■ 
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Buckets with Slater-determinants 

N=0 N=1 N=2 N=3 



step J 



copy 



step J+1 




FIG. 1: Visualization of the recursion step J ^ J + 1: for a considered bucket N copy all existing 
states, add new states by applying Cj to all existing states from previous bucket (A^ — 1) and 
truncate with the help of the measure /i. Filled circles represent single-particle states and N 
denotes the particle number in the bucket. (The shown truncation is just an example.) 

— "Expansion" 

Let Sn = Sj^N U {civ I V G Sj^n-i}- 

— "Truncation" 

If card(5'Ar) < Dn: Sj+i^n = Sn- 
Else: choose Sj+i^n C Sn such that 
card(S'j+i,Ar) = Dn 

and Vf G Sj^i^n '■ G Sn\Sj-,^i^n '■ l^{v) < fi{w). 
Here, |0) denotes the Fock space vacuum state, card(S') is the number of elements in set S, 
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and D]^{ (with N — 0, ...,Nmax) are given integers which determine the maximum number 
of states in Sj^n (for all J) . Note that the maximum number of Slater-determinants for N 
particles reads as 

TV I 

D (N) = (2) 

Hence, we must have Dn < DmaxiN). Furthermore, // is a suitable real "measure" for 
the selection of Slater-determinants. For a normalized Slater-determinant we assume a 
measure of the form //(?;) = {v\M\v) with a (w.o.l.g. diagonal) many-body operator M. The 
latter can be expanded in a series of n-particle product terms 

M = Y. M^l^N,, + Y: MfXNnNn + • • • , (3) 

31 jl<j2 

with Nj = CjCj for single-particle state j an real coefficients M^"'\ which represent conditional 
n-particle selection weights. In the simplest case, one can choose the diagonal elements of 
the Hamiltonian, that is 



Mj^) = (4) 

(2) 



^jk = ^jkkj - Vjkjk, 

M(">2) = 0. 



(Note that this choice might not provide a unique set if some Slater-determinants are degen- 
erate.) Optionally, the final truncation condition can be made even more selective, choosing 
for example only such Slater-determinants which have energies that are within certain energy 
intervals (such as "bands"). From this construction it is obvious that the set Sj^^ contains 
only Slater-determinants with a total number of particles. As a physical interpretation, 
with each recursion step J ^ J + 1 the set Sj+i^n consists of the lowest Dn states from 
the combination of Sj^n and all states of Sj^n-i with an extra particle created in the single- 
particle state j — J. To some extent, the described procedure resembles the NRG recursion 
scheme, however, lacking the diagonalization steps at this stage (see next paragraph). 

If we are solely interested in calculating expectation values of observables for a fixed 
particle number A^o, it is sufficient to calculate Sj^n for N = 0, ...,No only and to consider 
the last set SNmax,No in the following . However, if n-point Green's functions (i.e., contour- 
ordered correlation functions oin/2 creation and annihilation operators at n points in time) 
have to be calculated as well, we typically need information about all Fock subspaces of 
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particle numbers Nq — Nq — n/2, No + n/2. Obviously, n is restricted to < n/2 < 
mm{No, N^ax ~ Nq). Since the single-particle subspace basis Bj^^^^ might be optimized 
only for states with A^o particles, it could become sub-optimal for A^'o ^ Nq. In general, 
we thus could consider an individual subspace basis B^'^^_^ for each Nq and an individual 

Prom the obtained Sn^^^,n, we can now define a finite Fock subspace 

Fq = span(5') (5) 
spanned by the many-body basis S of "relevant" Slater-determinants 

U S^m..,N (6) 

N=No-^ 

and with dimension 



dim(Fo)= card(5^_,^)< ^ D^. (7) 

Note that the states of Si\imax,N for different N are orthogonal to each other since they 
obviously have different particle numbers. By construction, the states within Sn^^^^n for 
a fixed N are orthonormahzed. Hence, ,5 is a many-body ON-basis of Fq. In turn, a 
corresponding finite restricted Hamiltonian can be defined as 

Hq^ H \fo, (8) 

which acts within Fq. As the final step, Hq is diagonalized numerically in a matrix representa- 
tion with respect to the basis S. Since H conserves the particle number, H is block-diagonal 
within individual span(5'jv^^^^iv)- With the help of the eigenvectors and eigenvalues of Hq 
one can calculate the time evolution of the isolated system and can construct stationary 
statistical operators for the calculation of expectation values. 

By considering only those Slater-determinants that are chosen as "relevant" by the de- 
scribed algorithm, two controlled approximations are made: First, the number of considered 
single-particle states N^ax is finite. Equivalently, the recursion step is performed only a 
finite number of times Nj^^x- Second, the maximum number oi Slater-determinants 
in each set Sj^n can be smaller than Dmax{N). In the limit Nmax oo and Dat — > oo 
{N — 0, A^o), the algorithm becomes exact. 
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The described algorithm typically becomes advantageous if the required maximum num- 
ber of single-particle states Nmax and the number of particles Nq provides a number 
Dmax{No) of Slater-determinants that lies beyond the limit for which the resulting problem 
can be fully diagonalized, in particular for strongly inhomogeneous systems with minimal 
symmetries. Typical application-relevant values are N^ax = 32 — 256, iVg = — 32, and 
Dn = 1024 — 16384. For such a regime, NRG-based approaches might become impractical 
since the storage requirements of matrices for the given N^ax and Dn grow beyond realistic 
limits. In principle, the described algorithm is suitable to provide a basis for many-body 
states up to any energy scale since all single-particle states are scanned systematically. It 
can be used in particular for many-body problems where a single set of Slater-determinants 
has to be determined which is relevant for a whole energy interval of many-body eigenstates 
(not only the groundstate or a particular excited state). An implementation on parallel 
computer architectures is possible. 



III. OPTIMIZATION 



Within this algorithm, two means of (a-posteriori) optimization exist: (i) Choice of the 
single-particle basis and in particular its order, (ii) Choice of the selection criterion /i of 
Slater-determinants. 

As for (i), one could solve an effective single-particle mean- field problem derived from 
H for a giv en p article number Nq and the given occupation conditions (e.g., Hartree-Fock, 
LDA-DFT HQ)- 

The resulting single-particle ON eigenbasis of this problem could be 
chosen for B. Alternatively, one can define an outer self-consistency loop for the described 
algorithm. Here, the many-body state (preparation) of the system is decribed by a suitable 
many-body statistical operator p[H] which is expressed in terms of the calculated basis S 
and the associated eigenvalues for a given boundary condition (such as maximum entropy 
for thermodynamical equilibrium, non-equilibrium injection conditions, or the selection of 
a certain excited many-body state). In turn, the single-particle eigenstates of the resulting 
(transposed) single-particle density-matrix p, which can be calculated via 

Pkj = - G%{t = 0) (9) 
with G<(t) = «Tr(pcI(0)c,(t)), 
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are used as a new single-particle basis Bnew for a repeated many-body diagonalization, 
employing S[Bnew]- Such an approach resembles the MCSCF formalism [3] (where these 
single-particle states are referred to as "natural orbitals"), which is used typically for the 
calculation of the many-body groundstate. The correct choice of the single-particle basis is 
essential for any finite expansion of many-body states in Slater-determinants to be efficient. 

As for (ii), the measure /i can also be chosen adaptively from the many-body diago- 
nalization. Here, the coefficients M*^") can be constructed from the probability weights of 
conditional n-particle projection amplitudes of the calculated many-body statistical oper- 
ator p. Implemented as a repeated selection-|-diagonalization, this defines a self-consistent 
optimization scheme for /i. 

In a numerical implementation, the calculation (or transformation) and storage of the 
Coulomb matrix elements Vjkim might become a challenge {[{Nmax — '^)Nmax/'^]'^ independent 
index combinations). Here, symmetries of V have to be employed in combination with 
possible reduction and approximation schemes. 

IV. SUMMARY 

We have discussed a recursive method to construct a subset of relevant Slater- 
determinants for the use in many-body diagonalization schemes in order to simulate many- 
body effects of realistic nanodevices with Coulomb interaction. The system is assumed to be 
finite. The described method becomes advantageous if a large number of single-particle basis 
states (typ. 256) is required where full diagonalization schemes become impractical. By use 
of a recursive method, a reduced number (typ. a few 1000) of relevant many-body basis 
states (Slater-determinants) is generated, systematically scanning all given single-particle 
basis states. 
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